Renormalized transport properties of randomly gapped 2D Dirac fermions 
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We investigate the scaling properties of the recently acquired fermionic non-linear cr-model which 
controls gapless diffusive modes in a two-dimensional disordered system of Dirac electrons beyond 
charge neutrality. The transport on large scales is governed by a renormalizable nonlocal field 
theory. For zero mean random gap, it is characterized by the absence of a dynamic gap generation 
and a scale invariant diffusion coefficient. The j3 function of the DC conductivity, computed for this 
model, is in perfect agreement with numerical results obtained previously. 
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Transport in systems whose band structure has a node 
structure (e.g. Dirac points), as it appears in graphene 
and on the surface of 3D topological insulators, has been 
the subject of intense research recently. The experimen- 
tal observation of transport in graphene is character- 
ized by a minimal conductivity at the charge neutrality 
(or Dirac) point and by a linearly increasing conductiv- 
ity with increasing (electron or hole) charge density^. 
Thus, experimentally it is easy to distinguish the very 
robust minimal conductivity which is contrasted by the 
disorder-dependent conductivity away from the Dirac 
point. The behavior at the Dirac point has also been 
predicted by held theory, showing that the minimal con- 
ductivity is quite independent (or very weakly depen- 
dent) on disorder. For instance, a nonlinear sigma model 
approximation^— as well as perturbation theory in terms 
of disorder strength clearly shows a very weak disorder 
dependence^. This has led to the claim that ballistic 
transport cannot be distinguished from diffusive trans- 
port at the Dirac point. However, transport properties 
away from the Dirac point are theoretically not easily ac- 
cessible. There are several attempts, based on a classical 
Boltzmann approach, which predict the experimentally 
observed linearly increasing conductivity as we go away 
from the Dirac point. However, it was only recently that 
a more general held-theoretical approach, based on the 
Kubo formalism, was suggested to describe the transport 
of 2D Dirac fermions within a unified theory^ 7 -, using a 
four-body Hamiltonian. 

In this paper we focus on disorder due to a random 
gap. This case is particularly interesting because it can 
lead to a metal-insulator transition when the average gap 
is equal to a critical valued. Starting from a nonlinear 
sigma model that controls the diffusive modes, we study 
the renormalization of the interaction of these modes as 
well as the renormalization of the diffusion coefficient and 
the conductivity at the Dirac point and away from it. 



I. MODEL 

Below we give a brief sketch of the field theoretical ap- 
proach to the conductivity. Here we are led by the repre- 
sentation given in Refs. The main quantity to be 



computed is the disorder averaged two-particle Green's 
function. The disorder potential v of the strength g is 
supposed to have zero mean (v r ) = and Gaussian cor- 
relator (v r v r ,) — gS rr ,. For random gap disorder, the 
disorder averaged two-particle Green's function reads 



K rr , = -(Tr n [G rr ,(ie)a 1 Gj, r (ie)a 1 ]) v 



m.m'n.n' 



(i) 



where Tr n is taken on the extended Dirac space, <p is 
a four component superfield <f> — (ipi + , tpi _, x 2 X2,+)> 
consisting of a complex if>\ ± and a Grassmann X2 ± field. 
In our notation a 1 2 3 are usual Pauli matrices and a the 
2x2 unity matrix. The field averaging is defined as 



<•■•>* = / m 



with the action 

S = ■ (H + ie)4>) + g( 



(2) 



(3) 



where <7 3 
matrix, and 



H = 



l A ® <t 3 , e = el 4 <E) (T , 1 4 the 4x4 unity 

\ 



fH + fj, 
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H T - (i 

\ H T +n J 



(4) 



with the chemical potential /i, the non-random Hamilto- 
nian H n — ia ■ V and the random Hamiltonian 



H = H n 



VCTo 



(5) 



Then H is invariant under the global symmetry transfor- 
mation 



H = e S He § , 
where S is given by the following matrix 
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(6) 



(7) 
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with two scalar fields ip 1 and ip 2 , which obey Grassmann 
statistics, i.e. = -<^, and ip^- = —tpjipi- 

We decouple the interaction term in Eq. ((3]) by a 
Hubbard-Stratonovich transformation. Integrating out 
supcrficlds cj> yields an action in terms of composite su- 
persymmetric Hubbard-Stratonovich fields Q 

S' = iTrg(Q) 2 + logdetg(ff + ie + 2Qa 3 ). (8) 

A non-trivial vacuum of this theory Q is found from the 
saddle-point condition and turns out to be degenerated 
with respect to the transformation S defined in Eq. (J7J: 



3 e S = Q 1 + Q 2 e 



-2S 



(9) 



where Q 1 (Q 2 ) commutes (anticommutes) with S, and 
vanishes under the graded trace TrgQ = 0. On the 
saddle-point manifold; the action represents a fermionic 
non-linear er-model 



5 / = logdetg(ff 



2Q 1 a 3 + 2Q 2 a 3 e 2i >). (10) 



The field Q 2 represents the order parameter for the spon- 
taneous breaking of the symmetry generated by S. Ex- 
panding Eq. (fTUj) up to second order in Q 2 as 



S' =S Q +S", 
and using the exact relation for Grassmann fields 
e 2S = 1 + 2S + 2S 2 

yields 

S" = 4Trg [G Q 2 a 3 S 2 
+2(G Q 2 a 3 S) 2 +2(G Q 2 a 3 S 2 ) 2 } , 

with matrix Green's function 

G ,rr> =diag{ff+>ff-.0-)ff+} rr , ) 

and 

9±,rr< = [ H + W)<ra\rr>> 



(11) 
(12) 

(13) 
(14) 
(15) 



where rj ~ exp[— ir/g] is the scattering rate^i. Eventu- 
ally, we rewrite Eq. (|13[) in terms of scalar Grassmann 
fields (fj and obtain a nonlocal fermionic theory 



S[<p] 



9V 
2 



3=1,2 rr' 



J i' 



-2r) 2 (-iy Y Tr 2 {S9s,rr'9s,r'r} 



(p jr ,(p jr ,ip jr ip jr 



s=± 



which describes the diffusion of Dirac electrons, 
/i < rj, the diffusion coefficient reads 



D 



2^ + °^ 



,(16) 



For 



(17) 



II. RENORMALIZATION GROUP ANALYSIS 

Below we investigate the scaling properties of ac- 
tion Eq. (fl6|) at large distances. For this purpose we 
expand the Fourier transform of the vertex function 

Tr 2 {sg s rr ig s r / r } to the leading order in momenta 

s=± 

of the fields. This changes action Eq. (fT6|) to 



3,1 

~ iX Yl (- i y 5 k- q ,p-t{p- t ) 2 ^ikv'j q <Pj P v'ju (is) 

j,kqpt 

with the shorthands ^2 q = J d 2 q/(2Tr) 2 , 5k, P = 

(27r) 2 (5(fc + p), e' = gr/e/2, and D' = gr/D/2. The in- 
teraction strength is defined as 



A 



2 [if] 6 2 n 

3tt (rj 2 + /i 2 ) 2 3tt rj 



(19) 



The zeroth order term in field momenta of the interaction 
part is zero due to its locality, while the first order term 
vanishes by the symmetry. The frequency e is supposed 
to be small and sent to zero in the DC limit. For this 
reason we do not distinguish between e' and e. 

In order to find the infrared behavior of action Eq. (fT5)l 
we follow the usual prescription of the Wilson RG trans- 
formation^: We decompose Grassmann fields into fast 
ipj and slow ip s modes. The idea is to integrate out 
fast modes and to obtain an action which mimics action 
Eq. flTBl but contains solely slow fields. To the second 
order in DC perturbation theory this action reads 



S[<Psl ~ S [ip B ]+S int [ip s ] + (S int [ip s ,<p f ]) f 



DC 



■v{Swk['P a ><Pf] S i n t[ ( Pa><Pf])f 



DC 



(20) 



with 



-i\ (-iys k _ q , p _ t ( P -t) 2 

j.kqpt 



<P83k<P8j q <Ps3P<Psjt, 



(21) 



and S int [(p s , (pf \ represents terms which contain both slow 
and fast fields. The averaging operator reads 



/. . . \ dc _ ; 

\ // ^ 7 DC 



V[<p f ] 



e -s? l>t>f] 



with the free DC action 

r-DCr.„ l _ r>/ 



S° c l<p f ] = D'^V<p jrf V<p' j 



jrf- 



(22) 



(23) 



Obviously, this construction guarantees (l)f C = 1 and 
defines the DC propagator (in Fourier representation) 



^3 6 q,-k 

D'q 2 ' 



(24) 




DC 




FIG. 1: Diagrams responsible for renormalization of e. 



FIG. 3: Diagrams responsible for renormalization of A. 




FIG. 2: Diagrams responsible for renormalization of/?'. 



The derivation of the renormalization group equations 
for the frequency e, 'diffusion' coefficient D' and inter- 
action strength A of the action Eq. (fT8|) is a challenging 
task. Diagrams which have to be evaluated arise by merg- 
ing vertices shown in Figs. [TJ [5] and [3] The renormaliza- 
tion of the energy e comes from one-loop diagrams which 
emerge by averaging vertices depicted in Fig. [1] over fast 
fields. As shown in Appendix |A] the evaluation of the 
diagrams yields a result that does not develop any diver- 
gences in the infrared, since the vertex is proportional 
to the squared loop momentum, and the propagator to 
inverse squared loop momentum: 



2tt D' 



(25) 



where Ao denotes an upper cutoff. This expression is 
not a renormalization group equation in the strict sense, 
since it does not contain the running cutoff parameter 
£ = logA /A. It represents a kind of a finite size effect 
which disappears in the continuous limit. 

The renormalization of the diffusion coefficient can be 
obtained by integrating out fast fields in vertices depicted 
in Fig. [5J The evaluation of the functional integral is 
presented in Appendix [B] and yields the following renor- 
malization of the free action: 



]~f(p)<Pi P <Pi f 



[WtPirWir + ™?iM » (26) 



where m and D' are expansion coefficients of zeroth and 
second order in momentum p of the vertex function 



7(p) = 4 



D' 3 



E 

kq 



(k+p) 2 [(k+p) 2 - (q+p) 2 
k 2 q 2 (k + q + p) 2 



(27) 



There is no linear term in this expansion, since 7(p) is 
symmetric with respect to the sign mirroring of p, i.e. 
7(— p) = 7(p). In the Appendices IB II and IB 21 is shown 
that both coefficients m and D' vanish. This is a very 
important result, since it guarantees the preservation of 
the gapless diffusive mode and the realty of the diffu- 
sion coefficient even for the complex interaction strength. 
Therefore, the only running parameter is the interaction 
strength A. Its renormalization is due to the one-loop 



diagrams which emerge after integrating out fast fields 
in Fig. [31 Lengthy and elaborate calculations presented 
in Appendix [C] lead to the remarkably simple renormal- 
ization group equation 



7T D' 2 



(28) 



for each fcrmionic channel. This equation is easily solved 
with the same starting value A, = 
u given in Eq. (|T^|) 



in both channels, 



.(-I) 3 ' 



,2«2 



T 2 D ,4 



nD' 2 



7r 2 L>' 4 



(29) 



Both real and imaginary parts of the interaction scale 
down to zero but not equally fast. At large scales, the 
imaginary part of A becomes dominant and therefore gen- 
erates a genuine real interaction. The RG flow for the 
Grassmann field (p 2 is depicted in Fig. |U Fig. [5] shows 
the RG landscape in the parametric space spanned by 
the real and imaginary part of the interaction A^. The 
RG trajectories represent a set of excentric circles with 
diameter u , each attracted to the Gaussian fixed point 

This can be seen best if we 
in Eq. d28j) - Then we get 



at sftAj = Oand SsXj = 
substitute A^ = + ivj 



d eUj = 2u j v j , 

d tVj = (-iy (u 2 -v 2 ). 



(30) 



The right hand side of this system of differential equa- 
tions represents indeed a parametrized circle. 

Independently from the choice of the initial value both 
real and imaginary parts of A 2 (analogously for A x ) be- 
come equal at the length obtained from the condition: 



V * =1, with 4= log |, 



irD' 2 * ~° I 

with / denoting the mean free path, which gives 



£ = I exp 



irD' 2 



I exp 



32 p 



for n < rj. (31) 



Here we used definitions of the bare diffusion coefficient 
Eq. (TTTl) and interaction strength Eq. (flU)) . At half filling, 
i.e. for fi = 0, this scale is infinite. 



4 




FIG. 4: (Color online) Renormalization of the interaction 
strength A2. Solid (red) line shows the real part and dashed 
(blue) line the imaginary part of A 2 . The cross-over scale is 
shown by the vertical (black) dotted line at £« = itD^Uq 1 . 




FIG. 5: (Color online) Renormalization group flow of the in- 
teraction Aj in the parametric space spanned by u = KA^ and 
v — S$Aj. The flow in the upper halfplane corresponds to the 
Grassmann field ifi2, that in the lower plane to tpx. 



III. SCALING PROPERTIES OF THE DC 
CONDUCTIVITY 



Our ultimate task is to determine the scaling behavior 
of the DC conductivity. For this we need to compute the 
corrections to the conductivity which arise due to the 
doping. The DC conductivity is either determined from 
the Einstein relation a oc pD (p is the density of states 
at the Fermi level) or calculated from the Kubo formula, 



Here, the functional integral should be performed over 
the full action Eq. (HF 



(34) 



with (1) = 1. To the leading order in A^, the two- 
particles Green's function is approximated as 



K(q) 



■E 



(35) 



where 5 int is given in Eq. (|18l) . The functional integration 
is to be performed over the free action only. As shown in 
Appendix [D] we eventually arrive at the following one- 
loop RG equation for the conductivity 



y 3=1,2 



(36) 



Further progress can be made if we exploit Eq. ([29)1: 



4cr n 



ngD' 3 



1 



7T 2 D' 4 



(37) 



Unpleasant constants can be eliminated by rescaling 
u —> ttD' 2 u and using definition of the diffusion co- 



efficient D 1 = g/An. This finally yields 



'1 



u 2 l 



(38) 



The integration of this equation is simple and we obtain 
the following asymptotic expression for the conductivity 



a(u £)=a + -f\og[l 



(39) 



At half filling, i.e. u = 0, the conductivity does not flow, 
i.e. it is scale invariant. At large scales, i.e. for L £, £ 
given in ()31l) . the conductivity grows bi- logarithmically 
as function of the sample size u(L) ~ <r loglogL/£. For 
this reason £ can be associated with an intermediate lo- 
calization scale. Due to the infrared asymptotic freedom 
of the underlying model, this result should be asymptot- 
ically correct in all loops. 



IV. DISCUSSION AND CONCLUSIONS 



" 2 ^ m 



(32) 



The scaling properties of the conductivity are usually 
given by the /3 function 

P(a) = ^loga, (40) 



where the two-particles Green's function takes contribu- ^ ' dl 

tions from both channels j = 1, 2 into account: 

with rescaled logarithmic length I — u £. The conduc- 
tivity from Eq. (|39)) generates the /3 function as depicted 
K(q) = - (fiqf'jp)- (33) * n M Its shape reveals a striking resemblance of 

9 jj-p the numerically determined j3 function of graphene at 



5 




Appendix A: Renormalization of the energy e 

For sake of simplicity we use below the following nota- 
tion for the interaction part of the action Eq. (TT51) : 



FIG. 6: (Color online) /3 function corresponding to the con- 
ductivity in Eq. ([39)) . 



the Dirac point with random scalar potential disorder 11 : 
It starts at the value of the universal minimal conduc- 
tivity; it is strictly positive; it reveals a distinct maxi- 
mum related to the length £; it does not have any fixed 
points besides a = a and a = oo. Finally, it does 
not depend on any quantities apart from the conduc- 
tivity a itself, in line with the one-parameter scaling hy- 
pothesis. However, one might wonder whether the above 
mentioned bi-logarithmic asymptotics of the conductiv- 
ity compares well with the predicted logarithmic growth 
at the Dirac point^ i 10 ' 11 . It is indeed not difficult to 
reproduce the /? function of Ref. by applying the 
Wilson RG transformation directly to the two-particle 
Green's function and exploiting the scaling properties of 
the disorder strength g. Instead, in our approach we 
keep g scale invariant. This assumption suits well for 
weak disorder, provided the sample size is much smaller 
than disorder generated intrinsic length ~ exp[l/<7} 12 , 
which corresponds to the common experimental situa- 
tion. Under these circumstances the scale invariance of 
the conductivity was demonstrated both numerically and 
analytically 6 ? 1 ?. On the other hand, in finite samples 
the conductivity grows logarithmically as function of the 
chemical potential: a(u Q ) ~ a log(/x/7?) in accordance 
with Ref. [3. 

In conclusion, we have presented a scaling analysis of 
the diffusion coefficient and the DC conductivity of doped 
graphene with a random gap. For this purpose we have 
used an alternative held theory and investigated its scal- 
ing properties. On this basis we derived an invariant dif- 
fusion coefficient and an astonishingly simple expressions 
for the scaling of the conductivity that reproduces the 
distinct shape of the /3 function of disordered graphene, 
found previously in numerical calculations^. 



5 int = I dld2d3d4 6(1 + 3 - 2 - 4) 

x(3-4)Vjl^2^3^4- 

The renormalization of e is due to diagrams which are 
obtained by contracting fast fields in Fig. [TJ Contracting 
is possible in the only way, i.e. the diagrams have one- 
fold degeneracy. To take Grassmann statistics correctly 
into account it is necessary to permute fast fields through, 
such that they form 'normal ordered' pairs (ff')'. 

2iX^2(-l) j j dld2d3dA (3 - 4) 2 6(1 + 3-2-4) 
3 

= 2iA^(-l) J ' / dld2d3d4 (3-4) 2 6(l + 3-2-4) 
3 J 

The fields to be contracted are marked with black dotes. 
Then, the contractions can be performed and we obtain: 

2i\^(-l) j / dld2dMA (3-4) 2 (5(l + 3-2-4) 
3 J 

foi.^ 2 .*(3 - 4 ) n ( 4 ) - Pji.Vi4.*(2 - 3)H(3)] , 
where H(q) — 1/D'q 2 . First contribution is zero because 

d3 (3 - 4) 2 <5(3 - 4) = 0. 

Second contribution is finite and cutoff dependent: 

-2zA J d3 n(3)£(-iy J dl (3 - l)*<p jU <f/ jU 
3 

« -2 ? A J d3 H(3)(3) 2 ^(-iy j dl Vjls v' jls . 

The renormalization factor then reads: 

-2zA J d3 n(3)(3) 2 = -^A( A 2 -A 2 ), 

which reduces for A — > to — i\K^/2TiD' . Lifting it 
into the exponent and absorbing into the action gives 
Eq. (EH). 
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The renormalization of the diffusion coefficient comes 
from the diagrams constructed from vertices depicted in 



6 



Fig[2j Every diagram is twice degenerated, i.e. the ana- 
lytical expression reads 

A{-i\) 2 Y^{-^Y +3 J dld2d3dA J d5d6d7d8 



(3 - 4) 2 (7 - 8) 2 5(1 + 3 - 2 - 4)5(5 + 7-6-8) 

X [<P s li V 'f2i <Pf3i V>'f4i Vf5j ¥>'s6j <P f7j ¥>'f8j 

We permute fields through and perform functional inte- 
grations in order to get 

4(«A) 2 / dld2d3d4 J dbdM7d% <p lsi <p' 6sj 

ij 

x (3 - 4) 2 (7 - 8) 2 5(I + 3 - 2 - 4)5(5 + 7-6-8) 

X [( <Pf5 3 f'f2i ) ( ff3i f'fSj ){<Pf7j<P'f4i) 

= 4{iX) 2 Y^{-l) i+j 5 ij 6 ij 5 ij f dld2d3dA J d5d6d7d8 



\2/ 



x (3 - 4) 2 (7 - 8) 2 5(I + 3 - 2 - 4)5(5 + 7-6-8) 
<p lsi <p' 6sj 11(2)11(4)11(8) [5(5 - 2)5(3 - 8)5(7 - 4) 
-5(7-2)5(3-8)5(5-4)] 

= 4 ( a ) 2 E / dl <Pu</u J d2d4 (2 - f ) 2 

i 

x [(2 - i) 2 - (4 - 1 ) 2 ]n(2)n(4)n(2 + 4-1) 

i,p 

The function j(p) is defined in Eqs. and ([2"7]). The 
way, how momenta in arguments of 5-functions are in- 
tegrated out, is not unique. Therefore we can shift inte- 
gration variables when necessary. 



1. Computation of the mass term 



First we evaluate the expression for the mass: 



m oc 
d 2 k 



d 2 k f d 2 q k 2 [k 2 



(2tt) 2 J (2tt) 2 k 2 q 2 (k + q) 2 
d 2 q 



1 



q 2 (k + q) 2 (k + qf 



(2tt) 2 J (2tt) 2 
Shifting in the second term q q — k we have 
d 2 k r d 2 q 1 if d 2 k 

where I = log A /A. First term is conveniently evaluated 
using Feynman parametrization: 

d 2 k f d 2 q k 2 f d 2 k 



2tt) 2 J (2tt) 2 q 2 (k + qf 



(27T) 



x / dx 
Jo 



d 2 q 



1 



(2tt) 2 [(1 - x)q 2 + x(k + qf 



(*)■ 



Performing shift q ^ q — xk we get 



d 2 fc 



(2^) 2 
= 4^ ./ (2tt) 2 





d 2 k 



d 2 q 



(2tt) 2 [ 9 2 + x(I-x)fc 2 
f 



k 2 I dx 

£Hx(l — x ) 



dx 



f 

4^ Jo 

f r 1 ^ 

27T /„_/ X 



f 



d 2 fc 



x + 1 -x\ J (2tt) 2 
d 2 k £ f d 2 k 
27 



(27T) 



i.e. the very same result. Therefore, the mass is zero, the 
diffusive Goldstone mode is preserved. 



2. Computation of the diffusion coefficient 
renormalization 

Next we evaluate the renormalization of the diffusion 
coefficient: 



D' 



(i\) 2 l f d 2 q 



dp 2 



D' 3 2 
(q+p) 4 



d 2 k 



(2tt) 2 J (2tt) 2 

(q+p) 2 (k+p) 2 

k 2 q 2 (k + q+p) 2 k 2 q 2 (k+p + q) 2 



p=0 



We start with first term. Using Feynman parametriza- 
tion we have 



1 = 



dp- 



d 2 q 



d 2 k 



(q + pf 



(2tt) 2 J (2tt) 2 k 2 q 2 (k + q + p) 2 

= 8P_ f d 2 q (q+p) 4 '■' 
dp 2 J (2tt) 2 g 2 
d 2 /c f 



p=0 



dx 



(2ir) 2 [(1 - x)k 2 + x(k + q + p) 2 } 2 
Next shift k -+ fc — 01(3 + p): 

9 2 /■ d 2 g (q+p) 4 ^ 



p=0 



(27T)* 



(7x 



1 



47T <9p 2 

47T dp 2 

f 

47T 



<9p 2 

d 2 fc 

(2tt) 2 [fc 2 + x(I -x)(g+p) 2 

d 2 q (q+p) 4 
(2tt) 2 



p=0 



d 2 <j (<7+p) 2 



(27T) 



d 2 g 2 



(9 + p) 2 x(I - x) 

dx ~TT — T 

p= a Jo X \ L ~ x ) 

2 



p=0 



dx 



1 



(2^) 2 g 2 y 1(1 - x) 



It is important to recognize that shifting of the integra- 
tion variables with respect to the external momentum p 
does not affect the final result. This is because the inte- 
gration over the momentum conserving 5-function is not 



7 



unique. Indeed, if we shift q — > q—p in the above integral 
we obtain the same result: 



dp 2 



d 2 q 



d 2 k 



1 



(2tt) 2 (q-p) 2 J (2tt) 2 k 2 (k + q) 2 



p=0 



d 2 q 



2ir dp 2 J (2tt) 2 (q-p) 2 



p=0 



where we skipped integration over the momentum k. Af- 
ter mirroring sign of q we have 



£ d 2 



d 2 q 



q- 



2tt dp 2 J (2tt) 2 {q + p) 2 
I f d 2 q \ (e p -q) 2 ! 



2tt J (2tt) 2 



= (*)• 



Here we have to shed some light on the structure of the 
first term: 



d 2 q [e p -qf 
(2tt) 2 q 4 



l ° qdq 1 



271 



da(e p ■ q) z 



qdq 1 
v ° qdq 1 
qdq 1 



2tt 



<iag cos (a — (p) 
q 2 

da— (1 + cos 2(a — (^)) 



(2^) 2 2q 2 J 
Thus we continue: 



da 



d 2 q 1 
(2tt) 2 2g 2 ' 



(*) 



2^ 





'4 2 " 


/ (2tt) 2 


2 9 



d 2 g 1 
2vr J (2tt) 2 ^ 



= 2 



i.e. the very same result as before. This knowledge can 
now be used for evaluating second term: 



II = 



d 2 



d 2 q 



d 2 k (q+p) 2 (k + p) 2 



p=0 



dp 2 J (2tt) 2 J (2tt) 2 k 2 q 2 (k+p + q) 2 
Here we are allowed to shift k — > k — p, which yields 
d 2 q f d 2 k k 2 2 (q+p) 2 



(2tt) 2 J (2tt) 2 q 2 (k + q) 2 dp 2 (k - p) 2 

d 2 q f d 2 k k 2 I" 2 
(2tt) 2 J (2tt) 2 q 2 (k + q) 2 [k 2 
(e P ■ kf _ 2_\ (e p ■ q)(e p ■ k) 
k 6 k 4 ) k 4 
d 2 q f d 2 k 2 



p=0 



d 2 q 



(2tt) 2 J (2tt) 2 q 2 (k + q) 2 
Xe p -kf 



d 2 k 



1 



(27r) 2 y (2^) 2 [k + qf 



k 4 



2 

fc 2 



+8 



d 2 q 

(27T) 2 



<i 2 fc (e p • q)(e p ■ k) 
(2tt) 2 q 2 k 2 (k + q) 2 ' 



First and second terms can be easily evaluated after shift- 
ing k k — q and g — > q — k, respectively. Both give the 
same contribution I calculated above. The last term is 
more cumbersome: 



d q f d 2 k (e p ■ q)(e p ■ k) 



(2tt) 2 J (2tt) 2 q 2 k 2 (k + q) 2 



d 2 g e p ■ q 
(2tt) 2 q 2 



x / dx 



d 2 k 



,„ j (2ir) 2 [(l-x)k 2 +x(k + q) 2 } 2 ' 
Again, we shift k — > k — xq, which gives 



d 2 q e p ■ q f 1 f d 2 k 8e p ■ (k — xq) 



(2tt) 2 q 2 J J (2tt) 2 [k 2 + x{l - x)q 2 } 
d 2 q (e p ■ q) 



2 /■! 



dx 



d 2 k 



8.T 



(2tt) 2 q 2 J Q J (2tt) 2 [k 2 + x{l - x)q 2 } 2 
d 2 q (e p ■ q) 2 1 f 1 dx 



(2tt) 2 q 4 At: J 1-x 

1 f d 2 q 1 dy_ = _J£ 2 

ttJ (2n) 2 q 2 J 1 y \2n 

where at some point we have substituted y = 1 — x. 
Summing over all contributions to II we obtain 

II = I + I - I = I. 

The renormalization of the diffusion coefficient to order 
A 2 is therefore zero: 

D' oc I - II = 0. 



Appendix C: Renormalization of the interaction 
strength 

The renormalization of the interaction is due to one- 
loop diagrams which arise after merging vertices depicted 
in Fig. [31 We consider both contributions separately. 



1. Contribution I 
DC 




It has a two-fold degeneration, i.e. there are two possi- 
bilities for contracting fast fields: 

1= HA) 2 ^(-l) i+J J dld2d3dA J dbd6d7d8 

ij 

(3 - 4) 2 (7 - 8) 2 <S(1 + 3 - 2 - 4)5(5 + 7-6-8) 

[P s liV'f2iP s 3iV'fte9f5jPUjVf7j < P' S 8j 
+VsliV'f2iVs3iVwVf5jVs6j<Pf7j l p's8j]- ( C1 ) 
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Next, we permute and contract Grassmann fields: 
1= {i\) 2 Y^{-^) i+i f dld2d2,M j dbd6d7d8 
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(3 - 4) 2 (7 - 8) 2 5(1 + 3 - 2 - 4)5(5 + 7-6-8) 

+ (- 1 ) 1+5 ^sli(^f7j^f2'i} ( Ps3i( l Pf5j l P'f4i) i P's6j l P's8j] 

= {i\) 2 Yj(-l) i+i &i 3 bji J dld2dUA J dhd&dldZ 

ij 

(3 - 4) 2 (7 - 8) 2 5(1 + 3 - 2 - 4)5(5 + 7-6-8) 
x [5(7 - 2)5(5 - 4) - 5(5 - 2)5(7 - 4)] 

^i^ s3 y s6 y S 8,n(2)n(4). 

Performing summations and integrations we arrive at 
I = (iA) 2 J dld3d6d8 5(1 + 3-6-8) 
^V.uV^m'Pm / d4(4-3) 2 

i 

x [(4 - 6) 2 - (4 - 8) 2 ]n(4)n(4 - 3 - 1). 

After reordering fields in first term and renaming vari- 
ables we finally obtain 

I = -2(a) 2 J dld2d3d4S(l + 3-2-4) 

VnVaVaV'u / d5 (5 - 3) 2 

i J 

x(5 - 2) 2 n(5)IT(5 - 3 - 1). 

Now we expand vertex function up to the second order 
in fields momenta: 



/ 



d5 (5 - 3) 2 (5 - 2) 2 n(5)II(5 - 3 - 1) 



d 2 q (q-k) 2 (q-p) 2 



1 



D' 2 J (2tt) 2 q 2 (q-k-tf 



+ ■ 



^ 'dkdp 
t 2 _cP_ 

2 dt 2 



kt- 



d 2 



■pt- 



d 2 



k 2 d 2 



p^cf_ 

dkdt ' t '"dpdt ' 2* 2 dp 2 
1_ f d 2 q (q-k) 2 (q-p) 2 

r 2 J (2tt 



D' 2 J (2tt) 2 q 2 (q-k-t) 2 
Order kp: The factor is zero, since 
d 2 f d 2 q (q-k) 2 (q-p) 2 



(C2) 



k,p,t=0 



dkdp J (2tt) 2 q 2 (q - k) 2 

d 2 f d 2 q (g-p) 2 
J (2^ 2 



k,p=0 

= 



k,p=0 



dkdp J (2tt) 2 
after differentiation with respect to k. 
• Order kt: The expansion factor is also zero, be- 



cause: 



d 2 



dkdt 



d 2 q (q - kfq 2 



(2tt) 2 q 2 (q-k-i) 2 

d_ I" d 2 q 2e t -(q-k) 
~~dk J (2tt) 2 (q-k) 2 



k,t=0 



k=0 



d 2 q 



(2k) 2 
Rewrite first term: 

f2TT 



(*) 



/ d(p(e t ■ q)(e k ■ q) 
Jo 

= q 2 / d</?cos(o; — (/?) cos(/3 — ip) 
Jo 

q 2 f 27T q 2 f 2n 

= — cos(a - P) J dp = (e k ■ e t )— j dp. 
Plugging this back into (*) we see that 



d 2 q 



-I 



( 2 7t) 2 

d 2 q 



, (e k -q)(e t -q) _ ^ e t ■ e k 



(2k) 2 



o e fc • e t _ o et ' e k 



= 0. 



• Order pt: The expansion factor reads: 
d 2 f d 2 q q 2 (q-p) 2 



dpdt 



J (27T 



) 2 q 2 (q-t) 2 



p,t=0 



d 2 q d 2 d 

-Aq-p) 



1 



/ 



(2tt) 2 dp" 
d 2 q (e t -q)(e p -q) 
(2tt) 2 q 4 



dt (q - t) 2 



;t=0 



d 2 q e 



t ^ P 



(27T) 



Thus, the order pt in expansion is 
2tt 

Order k 2 : The corresponding factor is zero: 



1 d 2 r d 2 q q 2 (q-k) 2 
2dk 2 J (2tt) 2 q 2 (q-k) 2 



k=0 



1 d 2 
2W 2 



J (27T 



q_ 

) 2 



0. 



Order t 2 : The factor in expansion i 



is 



d 2 q q 2 



1 d 2 



2 dt 2 J (2tt) 2 (q - i) 2 



t=o 



Here, it is possible to shift q — > q + t and we obtain 



\_cP_ 

2dt 2 



d 2 q (q + t) 2 



(2tt) 2 q 2 



t=0 



2tt' 



• Order p 2 : The corresponding expansion coefficient 
reads 

a2 r ^ n „2f„ m \2 



1 d 2 r d A q 

2dp^ J (2^p 



q q 2 (q-p) 2 



p=0 



2tt' 
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Concluding, the leading order momentum expansion of 
the vertex function reads 

(t 2 -2t-p + p 2 ) £ (2-1) 2 I 



Eq.(G2 



D' 



2?r D' 2 2tt' 



and the contribution to the renormalization of the inter- 
action strength 

I ~ J dld2d3d4 6(1 + 3 - 2 - 4) 

x(2-l)Va^3^4- (C3) 

2. Contribution II 




Each diagram has the degeneracy one, and we may write 
corresponding expressions as 



II = 2(-i\) 2 J2(-iy +j J dld2d3d4 J d5d6d7d8 

ij 

(3 - 4) 2 (7 - 8) 2 6(1 + 3 - 2 - 4)5(5 + 7-6-8) 

+VsliV's2i <Pf3i V'fiiVsS] V'sdj V? f7j <Pf8j 
+^sli^'f2i^fm^Ui^s5j^m ( Pf7j ( Ps8j] ■ ( C4 ) 

After permuting fields, performing summations and in- 
tegrations, and renaming variables we obtain following 
expressions: 

Ha = 2(iA) 2 Y I dld2d'3d4 (3 - 4) 2 <5(1 + 3-2-4) 

xVutiiWu Jd5(5- l) 2 n(5)n(5 + 2-1), 
lib = -2(iA) 2 Y J dld2didA (3 - 4) 4 <5(1 + 3-2-4) 

i 

XtpiiVkVaitii J d5 11(5)11(5 + 3-4), 
lie = 2(iA) 2 Y J dld2d3d4 (3 - 4) 2 5(1 + 3-2-4) 

i 

xVutiiW'u J d5 (5 - 4) 2 n(5)n(5 - 2 - 4), 
lid = -2(iA) 2 Y J dld2d3dA 5(1 + 3-2-4) 

i 

xtPutiiVvtii Jd5(5~ 4) 2 (5 - l) 2 n(5)n(5 + 2-1). 



While evaluating contributions Ha and lie, it suffices to 
take only the most divergent part from the integral over 
loop momentum 5 into account. This yields 

IIa = IIc«2(zA) 2 !^g 2 n 2 (g) 

x Y J dld2d3d4 (3 - 4) 2 5(1 + 3 - 2 - 4)^^^ 



(zXf 



7T D' 



wYll dld2dsd4 5(1 + 3-2-4) 

i 

x(3~4) 2 ip li (p' 2l (p 3i (p' 4l . 



For contribution lb, it is necessary to perform a full inte- 
gration over the loop momentum with help of the Feyn- 
man parametrization: 



f d 2 q , . . . 1 f d 2 q 1 



ty 2 



_£ 1_ 

2^D'H 2 ' 



where we replace momentum (3 — 4) — > t. Hence, contri- 
bution lib reads 

lib = --^r Y J dld2d3d4: 8(1 + 3 - 2 ~ 4) 

i 

X(3 - 4) 2 tp u ip' 2l (p 3l tp' 4l = -Ha. 

The evaluation of the contribution lid goes analogously 
to the evaluation of the contribution I with the result 

lid = -2^f — V f dld2d3d4 6(1 + 3 - 2 - 4) 
D' 2 2n^J y ' 

x(2 + 4) 2 (p u ip 2i tp 3l tp 4i . 

An apparent problem with the momentum dependence 
of the vertex function can be cured if we remember that 
in DC limit ip and ip' are not independent but ip_ p = ip' . 
Using this property we can mirror momenta 2 — > —2 and 
3 — > —3 and obtain 



lid = 2 



(a) 2 i 



\ ^ 



J dld2d3d4 5(1 + 2-3-4) 



D' 2 2n^ 

x(2-A) 2 tp u <p' 3i ip 2i ip' 4i . 

The sing change is due to permuting Grassmann vari- 
ables. After renaming variables we get the topological 
structure of initial interaction term. Summing up all con- 
tributions gives 

I (iX) 2 

I + Ha + lib + lie + lid = -MrX 



7T D 



dld2d3dA (3 - 4) 2 5(1 + 3 - 2 - 4)^^^ 
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Therefore, the RG equation for the interaction strength 
acquires the form 



V 7T D 



(C5) 



It is therefore convenient to distinguish between interac- 
tions in each channel: 



i\ 1 = iX, + (-1) 



1V^ 



7T D' 



In continuous limit this gives Eq. 



Appendix D: Scaling of the DC conductivity to 
one— loop order 

The DC conductivity is calculated from Kubo formula 
Eq. (I3"2"j) with the two-particles Green's function defined 
in Eqs. (153")) and The action is slightly changed 

using the acquired knowledge, as 



i.e. the usual universal DC conductivity of Dirac electron 
gas. The evaluation of the second term is more cumber- 
some. Respecting all possible (four in total) contraction 
combinations yields: 

*i(9)~-E / MiAM) 
a J 

ija 

J dld2d3d4 (3 - 4) 2 5(1 + 3 - 2 - 4) 

+Viq<P'jpV a l V' a 2'Paa<Pa4 + <Piq<Pjp a l t Pa2'P a 3^' a 4 



a 4 
'A I ' 



ija 



4) 



with 



5 oM = E / dld2 bA 1 - 2 ) + D 'v 2 W, 

a a i n " 



and 



/ (2^ / dld2dMA ( 3 ~ 4) 2 (5(1 + 3 - 2 
[-6(2 - q)5(l - p)S(3 - 4)K Q (q)K Q (p)K (3) 
+5(4 - q)S(l - p)5(3 - 2)K Q (q)K Q (p)K Q (3) 
-5(4 - q )5(3 - p)6(l - 2)K Q (q)K Q (p)K Q (2) 
+5(q - 2)5(3 - p)5(l - 4)K (q)K (p)K (4)} . 

While contributions from the first and third terms vanish, 
both other give equal finite contributions: 



Si„tM = E (-i)^ / dld2d3d4 5(1 + 3-2-4) K ^ ~ 2% E {D > q 2 + / e) 2 J (2?r)2 W + le 



d 2 k e 



(D2) 



i=i,2 



x (3-4) 2 (p J - 1 ^- 2 ^3^ 4 . 
We evaluate Eq. (|35[) denoting 

-^1(9) ~ E^i9^"P 5 int) ' 



The integral 



d 2 /c 



(2tt) 2 D'k 2 



where 



behaves well in DC limit (e — > 0) as it does not develop 
any IR divergences. The corresponding contribution dis- 
appears in continuous limit and does not affect the con- 
ductivity at large scales. Therefore, the main contribu- 
tion arises from the following expression: 



<->° 



Zn 



DU... e -5 M. 



KM 



2iq 2 
(ief 



E 

.7=1,2 



d 2 k [(-l) j Xj 



(2tt) 2 D'k 2 + 



Here, the proportionality factor is 1/g. The zeroth order 
contribution to the conductivity is calculated as usual: 

^)~E/(02(^)° 



In DC limit, the integral diverges logarithmically. The 
corresponding conductivity correction reads 



a 1 = 2e 



2 ^ KM 



2 ^ E [(-^] ^ 



E 



d 2 p <5„ (2tt) 2 5(q - p) 



(2n) 



D'q 2 + it 



D'q 2 



and further with D' — g/4n: 
e 2 d 2 



i=o * 3=1,2 

and the full expression for the renormalized conductivity: 

(D3) 



y 3=1,2 



a = 2- 



g dq 2 D'q 2 



2D' _ 1 

9 =o 9 v 



(Dl) Continuous limit of this expression yields Eq. (|3l))) . 
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